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ABSTRACT 

We present a semi-analytical model of star formation which explains simultaneously the 
observed UV luminosity function of high redshift Lyman break galaxies (LBGs) and lumi- 
nosity functions of Lyman-a emitters. We consider both models that use the Press-Schechter 
(PS) and Sheth-Tormen (ST) halo mass functions to calculate the abundances of dark matter 
halos. The Lyman-a luminosity functions at z < 4 are well reproduced with only < 10% of 
the LBGs emitting Lyman-a lines with rest equivalent width greater than the limiting equiv- 
alent width of the narrow band surveys. However, the observed luminosity function at 2; > 5 
can be reproduced only when we assume that nearly all LBGs are Lyman-a emitters. Thus it 
appears that 4 < z < 5 marks the epoch when a clear change occurs in the physical properties 
of the high redshift galaxies. As Lyman-a escape depends on dust and gas kinematics of the 
inter stellar medium (ISM), this could mean that on an average the ISM at z > 5 could be less 
dusty, more clumpy and having more complex velocity field. All of these will enable easier 
escape of the Lyman-a photons. At z > 5 the observed Lyman-a luminosity function are well 
reproduced with the evolution in the halo mass function along with very minor evolution in 
the physical properties of high redshift galaxies. In particular, upto z = 6.5, we do not see the 
effect of evolving inter galactic medium (IGM) opacity on the Lyman-a escape from these 
galaxies. 
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1 INTRODUCTION 

Determining the star formation history of the high redshift universe 
is one of the major goals of ongoing observations. Available obser- 
vational data mainly consists of UV luminosity functions (LFs) of 
high redshift Lyman break galaxies (LBGs) which can in turn give 
the star formation rate density of the universe. The galaxies have 
been identified even up to redshift z ~ 10 using so called photo- 
metric 'drop-out' technique (Bouwens et al. 2004, Hopkins & Bea- 
com 2006, Richard et al. 2006). However, very good constraints are 
available only up to z ~ 7 (Bouwens et al 2008). 

In addition to the 'drop-out' techniques, narrow band searches 
for high redshift galaxies emitting a strong Lyman-a line are suc- 
cessful in detecting galaxies at 3 < 2 < 6 (Cowie & Hu 1998, Hu 
et al. 1998, Rhoads et al. 2000, Taniguchi et al. 2005, Shimashaku 
et al. 2006, Kashikawa et al. 2006, Murayama et al. 2007, Gron- 
wall et al. 2007, Dawson et al 2007, Ota et al. 2008, Ouchi et al. 
2008). Unlike the drop-out technique used in detecting the LBGs, 
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the searches for Lyman-a emitters are not biased by UV luminos- 
ity. However, the detectability depends on the Lyman-a emissivity 
and radiative transport. Thus these two techniques pick up galax- 
ies with different types of selection biases. Availability of the UV 
luminosity functions of Lyman-a selected galaxies allows us to un- 
derstand these biases and provides joint constraints on models of 
galaxy formation at z > 3. 

The star formation rate is a key quantity for both UV as well as 
Lyman-a emission from a galaxy. Hence, it is interesting to obtain 
a semi-analytical model of star formation for these high redshift 
galaxies that can explain both these sets of observations. In our 
previous work by Samui, Srianand & Subramanian (2007) (here- 
after Paper I) we have built a semi-analytic model of star formation 
taking account of several feedback processes in order to explain 
the observed UV luminosity functions of LBGs at 3 < 2: < 10. 
By fitting the observed data we put constraints on the nature of the 
star formation in this redshift range. In Samui, Subramanian & Sri- 
anand (2009) (Paper II) we studied the effect of assumed form of 
the halo mass function on the results of semi-analytical galaxy for- 
mation models, in detail. As a continuation of these works, here 
we compute the luminosity function of Lyman-a emitters (LAEs) 
using the same star formation model and compare it with the three 
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sets of available observations, which are the high redshift UV lumi- 
nosity functions of LBGs and UV & Lyman-a luminosity functions 
of LAEs. Previous semi-analytical works on high redshift luminos- 
ity functions of Lyman-a emitters (i.e. Haiman & Spaans 1993, 
Thommes & Meisenheimer 2005, Le Delliou et al. 2005, 2006, 
Kobayashi et al 2007, Mao et al 2007, Dijkstra et al 2007, Stark 
et al. 2007) have considered a more limited set of currently avail- 
able observations. Our present work using this more extensive set, 
i.e. UV and Lyman-a luminosity functions of LAEs and UV lumi- 
nosity functions of LBGs in the redshift range 3 < z < 6.5, allows 
us to constrain the physical properties of LAEs and their redshift 
evolution. We use the cosmological parameters consistent with the 
recent WMAP data (Dunkley et al. 2008) = 1, = 0.26, 
= 0.74, fit = 0.044, h = 0.71, ag = 0.80 and Us = 0.96). 



2 SEMI-ANALYTICAL MODELS 

In order to compute the luminosity function of high redshift galax- 
ies one needs to model both the star formation in an individual 
galaxy and the abundance of dark matter halos in which the galax- 
ies form. We compute the abundance and formation rate of dark 
matter halos as a function of redshift in the framework of Lambda 
cold dark matter (LCDM) cosmology. For this purpose we con- 
sider two halo mass functions, the analytically motivated Press- 
Schechter (PS) halo mass function (Press & Schechter, 1974) and 
the Sheth-Tormen (ST) halo mass function (Sheth & Tormen 1999), 
which gives a better fit to numerical galaxy formation simulations. 
For the PS halo mass function, we use the formalism of Sasaki 
(1994) to calculate the net formation rate of halos. However, the 
Sasaki formalism is not easily genaralisable to the other form of 
mass functions. Hence for the ST halo mass function, we simply 
take recourse to its derivative to calculate the net formation rate of 
dark matter halos (also see Paper II). 

The star formation rate of an individual galaxy of dark matter 
mass M is assumed to be (Chiu & Ostriker 2000), 



Msf(M,z,Zc) = /* 



fib 



M 



t{z) - t{Zc) 



X exp 



t{z) - t{Zc) 



(1) 



^ ^dyn {^c, 

where, the amount and duration of the star formation is determined 
by the values of /* and k respectively. We can fix these two pa- 
rameters by fitting the observed UV luminosity functions of high 
redshift LBGs (see Paper I and II for details). Further, t{z) is the 
age of the universe; thus t{z) — t{zc) gives the age of the galaxy at 
z that has formed at an earlier epoch Zc, and tdyn is the dynamical 
time at that epoch. The star formation rate is converted to Ivmiinos- 
ity at 1500 A assuming an initial mass function (IMF) of the stars 
formed (see Eq. (6)-(8) of Paper I). The observed luminosity is less 
by a factor, rj, than the actual luminosity because of the dust red- 
dening inside the galaxy. In principle, the value of 7] depends on 
the wavelength and the functional form is governed by the nature 
of the dust grains. As in Paper I, we calculate the reionization his- 
tory of the universe and the radiative feedback of the meta-galactic 
background UV radiation on the star formation in a self-consistent 
manner for each model (also see Thoul & Weinberg 1996, Bromm 
& Loeb 2002, Benson et al. 2002; Dijkstra et al. 2004). We assume 
a steep cut off of star formation in halos with mass 
which is attributed to AGN feedback (Bower et al. 2005; Best et al. 
2006). 

We compute the Lyman-a luminosity of a star forming galaxy 



assuming case-B recombination. In this case, two Lyman-a pho- 
tons are produced out of three hydrogen ionizing photons (Oster- 
brock, 1989) that are confined within the interstellar medium of the 
galaxy. Hence the Lyman-a luminosity produced in any star form- 
ing region is related to its star formation rate by. 



Llvcc = 0.68/ii^„(l - Uc)N-fMsF. 



(2) 



Here, hva = 10.2 eV and /esc =0.1 are the energy of a Lyman-a 
photon and the escape fraction of UV ionizing photons respectively. 
Further, A^^. is the rate of ionizing photon production per unit solar 
mass of star formation. This mainly depends on the initial mass 
function of the stars and also on the metallicity. Values of number 
of ionizing photons per baryon of star formation for different IMFs 
and different metallicities can be found in Table 1 of Paper I. The 
observed Lyman-a luminosity is given by 

^LyoL — /esc ^Lyot- (3) 

Here, /<^c" is the escape probability of the Lyman-a photons. This 
is decided by the dust optical depth, velocity field of the ISM in 
the galaxies and the L3Tnan-a optical depth due to ambient inter- 
galactic medium around the galaxies. As Lyman-a is a resonant 
transition we expect the effective dust optical depth for Lyman- 
a in the ISM to be much larger than that for the UV continuum 
photons (i.e /^c" < I/'?)- However, if Lyman-a emission comes 
from some outflows in the star forming region (Malhotra & Rhoads 
2002, Dijkstra et al. 2007a, Verhamme et al. 2008) or through in- 
homogeneous ISM (Neufeld 1991, Hansen & Oh 2006, Finkelstein 
et al. 2008, 2009) then there may not be any correlation between 77 
and/^J". 

The escape fraction of Lyman-a also depends on the optical 
depth of the IGM in the immediate neighbourhood of the galaxy, in 
particular the proximate region that is affected by excess ionization 
by the galaxy itself. Thus the redshift evolution of /^c" can be an 
useful probe of the reionization history of the universe (Malhotra 
& Rhoads 2004, Stem et al. 2005, Haiman and Cen 2005, Dijk- 
stra, Wyithe & Haiman, 2007) and/or the redshift evolution of dust 
abundance (Mao et al. 2007), velocity field and gas clumping factor 
in galaxies. 

It may be possible that all the LBGs do not have a detectable 
Lyman-a emission. The spectroscopic observations of LBGs by 
Shapley et al. (2003) show that only 25% of the LBGs at z ~ 3 have 
Lyman-a emission with rest equivalent width Wq > 20 A(also see 
Steidel et al. 2000). Also the observations of UV luminosity func- 
tion of Lyman-a emitters show similar results. Hence we consider 
that only a fraction Gf of the entire galaxy population will be de- 
tected as Lyman-a emitters in surveys as they are usually sensitive 
to galaxies having Lyman-a equivalent widths above certain limit- 
ing value. 

The rest frame equivalent width of the Lyman-a emission is 

given by 



Wo = Ll^y^/iLcont/ri) = f^c°'LLyc/iLcont/ri) 



(4) 



where Leant is the continuum luminosity per unit wavelength 
near 1215 A. We obtained this from the stellar synthesis code 
'Starburst99'^ (Leitherer et al. 1999). For our continuous mode of 
star formation we use the same prescription as in Paper I for the UV 
continuum flux, to calculate the 1215 A continuum flux. We tab- 
ulate the continuum luminosities at different wavelengths and rest 
frame equivalent width of Lyman-a emission in Table I for various 
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Table 1. Flux at various wavelengths as predicted by 'Starburst99' for a continuous mode of star formation at a rate of 1 Mq/jt. The quoted values ai'e at the 
time t = 1.4 X 10® yrs after the star formation began. We also show the equivalent width calculated at t = 10^ yrs. 



Model IMF Metal /J^g /j^g^g ratio HI-UV^ /j'jig EW* (A) 

/905//1505 1.4 X 10»yrs 1.0 X lO'yrs 



model 1 


1 


100 


0.050 


40.12 


40.37 


0.57 


53.13 


40.44 


49.4 


52.4 


model2 


1 


100 


0.040 


40.10 


40.32 


0.60 


53.23 


40.36 


74.3 


81.2 


models 


1 


100 


0.020 


40.15 


40.42 


0.53 


53.33 


40.45 


75.8 


87.0 


model4 


1 


100 


0.008 


40.16 


40.46 


0.50 


53.45 


40.45 


100.3 


122.9 


models 


1 


100 


0.004 


40.18 


40.48 


0.51 


53.50 


40.46 


110.6 


137.9 


model6 


1 


100 


0.001 


40.21 


40.48 


0.53 


53.57 


40.44 


135.3 


169.6 


model? 


1 


100 


0.0004 


40.24 


40.50 


0.54 


53.61 


40.48 


133.4 


173.6 


models 


0.1 


100 


0.0004 


39.83 


40.10 


0.54 


53.20 


40.07 


139.4 


173.3 


model9 


5 


100 


0.0004 


40.58 


40.75 


0.66 


53.94 


40.80 


135.5 


177.6 


model 10 


10 


100 


0.0004 


40.75 


40.79 


0.91 


54.11 


40.92 


157.6 


181.8 


model 1 1 


20 


100 


0.0004 


40.89 


40.78 


1.29 


54.29 


40.95 


216.7 


216.7 


model 12 


40 


100 


0.0004 


40.95 


40.77 


1.54 


54.39 


40.93 


288.8 


288.8 


^ log of flux (erg s ^ 


at 905 A. 
















^ log of flux (erg s^^ 


A-i)at 1505 A. 

















t t Note that 'Starburst99' gives flux at 905 A and 1505 A. 
^ log of no. of H I ionizing photons per sec. 
4 log of flux (erg s"! k'^) at 1215 A. 

* Equivalent widths are calculated taking f^s^ = 1 and fesc = 0.1 and rj = 1. 
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Figure 1. The time evolution of intrinsic Lyman-cj equivalent width of 
a galaxy as predicted by our models. Solid and dashed lines are for the 
Salpeter IMF in the mass ranges 1 — 100 Mq and 10 — 100 Mq respec- 
tively (i.e. model7 and modellO of Table 1). We have assumed k = 1. The 
dynamical time (tdyn) depends on redshift of collapse (Zc). For example, 
at Zc = 10, 5 and 3, idyn = 8.8 X lO'^ , 2.2 X 10** and 4 X 10** yrs 
respectively. We also show the actual time that corresponds to Zc = 10 on 
the top axis. 



physical parameter related to the nature of the star formation for a 
continuous constant star formation model as obtained from 'Star- 
burst99' at the time t = 1.4x10* yrs after the star formation began. 
We also quote the equivalent width calculated at i = 10^ yrs. Note 
that for a constant continuous star formation model, the number of 
high mass star remain constant after typical life time of OB stars 
that dominate in Lhya- Hence, the equivalent width decreases with 
time as contribution from the low mass stars to L^ont continuously 
adds up. 



In our model the Lyman-a equivalent width of a galaxy is in- 
dependent of its mass as both Lyman continuum as well as line 
flux would scale with mass. It depends on the value of ntdyn and 
most importantly the values of rj and f^Ji" ■ It also depends on the 
value of /esc, the escape fraction of the ionizing photons. In Fig. 1 
we show the time evolution of the intrinsic rest frame equivalent 
width of a galaxy as predicted by our models. The observed equiv- 
alent width would be scaled by a factor rj/^Jf.". We also assume 
a Salpeter IMF in the mass range 1 — 100 A'Iq (solid line) and 
10 — 100 AIq (dashed line). As can be seen from Table 1 as well 
as from Fig. 1, the intrinsic equivalent width depends on the as- 
sumed IMF. Note that through out this work we will use k — 1 and 
Salpeter IMF with the mass range 1 — 100 A^q with metallicity 
0.0004 (i.e. model? of Table 1). 

Note that we have mainly three sets of observations that can 
be used to constrain our model parameters. These are (i) UV lu- 
minosity function of LBGs, (ii) Lyman-a luminosity function of 
Lyman-Q emitters and (iii) UV luminosity function of Lyman-a 
emitters. Along with these we have the information about the equiv- 
alent width distribution of the LAEs. The first set of observations 
can be used to constrain f*/r) combination. The second set can be 
used to constrain /, f^sc"' and the last one can be used to obtain G / . 
Then we will be able to calculate the mean Wo. The spread in the 
equivalent width of the detected galaxies will come in two ways: (i) 
distribution in and /^^c" and (ii) the spread in their ages. Since in 
our model we assume only the average value for both ri and /<ic°, 
we will have distribution in Wo only coming from the spread in the 
ages of detected galaxies. We show this distribution in the follow- 
ing section while discussing our results. 



3 LUMINOSITY FUNCTIONS 

In Fig. 2 we compare our model predictions for both UV and 
Lyman-a luminosity functions with the observed data points. For 
each redshift bin we have used the most recent measurement of the 
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Figure 2. Top panels: The predicted UV luminosity functions of LBGs (thin green lines) and LAEs (thick blue lines) at ^ = 3,4 and 6 along with the observed 
data points. The observationally determined UV luminosity functions of LBGs are taken from Reddy & Steidel (2008) (open red circle at 2 = 3) and Bouwens 
et al (2007) (open red circles at ^ = 4 and z = 6). The black filled and open triangles are observed data points of UV LPs of LAEs from Ouchi et al. (2008) 
for the reliable and less reliable points respectively. We show our model predictions for both Press-Schechter (dashed lines) and Sheth-Tormen (solid lines) 
halo mass functions. Bottom panels: The predicted Lyman-Q LF of LAEs at 2 = 3.1, 3.7 and 5.7. The observed data points (red filled triangles) are from 
Ouchi et al. (2008). The solid lines are for the models with the Sheth-Tormen mass function where as dashed lines are for the Press-Schechter mass function. 



Table 2. Comparison of model predictions between the PS and ST mass functions. 



PS Mass function ST Mass function 







Gf 


f fLya% 
J^Jesc 


EW* (A) 




Gf 


f fLya% 
J * J esc 


EW* (A) 


3.1 


0.044(3.90) 


0.07 


0.059 ±0.011 (0.95) 


179 


0.055(0.97) 


0.07 


0.076 ±0.011 (2.49) 


183 


3.7 


0.046(2.32) 


0.10 


0.051 ±0.014(0.68) 


148 


0.042(1.09) 


0.10 


0.050 ±0.015 (1.08) 


159 


5.7 


0.081(1.19) 


1.00 


0.044 ±0.017 (0.91) 


72 


0.050(0.63) 


1.00 


0.028 ±0.021 (0.42) 


75 


6.5 




1.00 


0.054 ± 0.012 (2.30) 






1.00 


0.031 ±0.015(2.32) 





t obtained using minimization and also coiTected for dust opacity at A = 1500 A ; 

the per degree of freedom are given in bracket (see Paper II for details), 
t values indicated inside the bracket are best fit x^ per degree of freedom. 
* the average equivalent width is calculated at t = 10** yrs. 



UV luminosity function of LBGs that covers a wide range in lumi- 
nosity. Below we provide details of observational data used in each 
redshift bins. Luminosity functions of LAEs are taken from Ouchi 
et al. (2008). The solid and dashed lines are our model predictions 
using ST and PS halo mass functions respectively. 

At a particular redshift, we first fit the observed UV lumi- 
nosity functions of LBGs by adjusting /,/r;. For this we use 
minimization technique (see Paper II for details). Then we fit the 
observed UV luminosity function of LAEs by changing G / and 
keeping same /* /t; obtained for the nearest available redshift. Note 
that, we did not try to get G / through minimization as there are 



only few data points in the observed luminosity function (also there 
are issues related to the completeness of the samples). Finally we 
match our model predictions with observed Lyman-Q luminosity 
function at the same redshift by adjusting ftfcsc"' and keeping Gf 
fixed. This is also done using minimization. In Table 2 we sum- 
marize the best fit parameters along with the x^ vales at different 
redshifts for models with both PS and ST mass functions. Below 
we describe these results for specific redshifts. 
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EW(A) 

Figure 3. Distribution of Lyman-o equivalent width as predicted by our 
model at 2 = 3.1,3.7 and 5.7. The thin lines are for the models that assume 
the PS halo mass function where as the thick lines are for the models with 
the ST mass function. Note that, here the spread in EW only comes from 
different ages of galaxies contributing to the luminosity function. In reality 



Lya 



more spread is expected from the spread in values of ri and /, 



3.1 Luminosity functions at z ~ 3 

At z '-^ 3 we have all the three observed luminosity functions and 
the data are quite well established, i.e. different groups have con- 
firmed the data by different methods. For 2; ~ 3 we use the ob- 
served UV luminosity function of LBGs given by Reddy & Steidel 
(2008) which covers the low luminosity end well. We show in the 
left most panels of Fig. 2 (panel (a) and (b)), both UV and Lyman-a 
luminosity functions at 2; ~ 3 as predicted by our model along with 
the observed data points. The best fit model parameters are given in 
Table 2. In panel (a) of Fig. 2, the set of thin curves in the top are 
the predictions of UV luminosity functions of LBGs. 

A good agreement with the observed UV luminosity function 
of LBGs is obtained for /./t? = 0.044 ± 0.001 and 0.055 ± 0.001 
for models with the PS and ST mass functions respectively. The 
corresponding reduced for these fits are 3.9 and 0.97. Thus the 
shape of the observed UV luminosity function of the LBGs is bet- 
ter reproduced by the model with the ST halo mass function. We fit 
the UV luminosity function of the Lyman-a emitters by multiply- 
ing the UV LF of LBGs with a fraction G/ . As mentioned earlier, 
we did not try to get G/ through minimization. The observed 
data points are well reproduced for Gf — 0.07 with same f,/r]. 
These curves are also shown in the figure by the thick blue lines 
(the bottom set of curves in panel (a)) for both PS and ST mass 
functions. The declining trend in the UV luminosity function of 
Lyman-Q emitters seen in the low luminosity end (open triangles) 
is mainly due to incompleteness. Apart from these points other data 
points do not require luminosity dependent Gf. This is consistent 
with our implicit assumption that G/ is independent of halo mass 
(or galaxy luminosity). 

The value of G/ is in agreement with the measurements of 
Shapley et al (2003) where they found that the fraction of LBGs 
having Lyman-a emission with equivalent width Wo > 60 A is 
~ 8% (see their Fig. 8). Note that the sample of LAEs of Ouchi et al 
(2008) at 2 = 3.1 has Wo > 60 A. Hence our results are consistent 
with both these observations. From the fact that G/ matches with 
the prediction from fig. 8 of Shapley et al. (2003) we can conclude 
that both the techniques of detecting z ^ 3 galaxies appear to pick 
a subset of the same parent population of galaxies. 



We now turn to Lyman-a LF of LAEs. We show this in 
panel (b) of Fig. 2 along with the observational data taken from 
Ouchi et al. (2008). To fit the Lyman-a luminosity function, only 
free parameter is as G/ has already been fixed by fitting 

the UV luminosity function of the LAEs. The best fit with the ob- 
servational data are obtained with /./^^r = 0.059 ± 0.011 and 
0.076 ± 0.011 for the PS and ST mass functions respectively. The 
corresponding best fit are 0.95 and 2.49. If we consider rj — 4.5 
as obtained by Reddy et al. (2006) then we have f^sl" = 0.29 and 
0.30 respectively for the models with the PS and ST mass func- 
tions. Taking rj = 4.5 also implies /, — 0.20 and /* = 0.25 for 
the models with the PS and ST mass functions respectively. 

We now calculate the average rest frame equivalent width of 
the Lyman-a emission of the star forming galaxies that are con- 
tributing to the luminosity function. Note that for given values of 
/*/j7 and f*ftsc°', the equivalent width is solely determined from 
the IMF we assume. This can be easily understood if we rewrite 
Eq. 4 as 



Wo = 



(5) 



Lcont{f*/ri) 

The ratio L^ya/ Leant depends on the IMF and the metallicity of 
the gas (see Table 1) and f*/ri and /*/ic" come from the fit. 
Therefore, fitting simultaneously the UV and Lyman-a luminos- 
ity functions uniquely specify the average equivalent width of the 
Lyman-a emission line. For the fit presented in panels (a) and (b) 
of Fig. 2 the average equivalent widths are 179 A and 183 A for 
the models with the PS and ST mass functions respectively. Note 
that 'rj' reflects extinction at A ~ 1500 A; the relative extinction at 
A = 1215 A will be higher than that at A = 1500 A. Therefore, the 
actual equivalent width will be higher depending upon the adopted 
extinction correction. The Lyman-a rest equivalent width distribu- 
tion predicted for the best fit model parameters is shown in Fig. 3. 
Note that the spread in rj and f^sc" around their best fitted val- 
ues will make this distribution spread over wider equivalent width 
range. Hence, one should not directly compare this histogram with 
observations although the mean value itself is relevant. In all re- 
sults presented here we use a lower mass cut off of 1 Mq in the as- 
sumed IMF. Increasing this to > 10 A/q to mimic a top-heavy IMF 
would increase the predicted equivalent width by a factor ~ 1.4 
(see Fig. 1). 

3.2 Luminosity functions at z ~ 4 

We show our model prediction as well as the observed data points 
at z ~ 4 in panel (c) and (d) of Fig. 2. The observed UV luminosity 
function of LBGs at 2 = 4 is taken from Bouwens et al. (2007). The 
luminosity function of Lyman-a emitters is at 2 = 3.7 and we com- 
pare this with UV luminosity function of LBGs at 2 = 4. We see 
from our model predictions that there is no significant change in the 
properties of the galaxies from 2 ~ 3 to 2 ~ 4. The UV luminosity 
function of LBGs can be well fitted with = 0.046±0.001 and 
0.042 ±0.001 with all other parameters being same as at 2 = 3 for 
the models using the PS and ST mass functions respectively (see 
Table 2). The corresponding reduced are 2.32 and 1.09. Thus, 
even for this redshift bin the model with the ST mass function pro- 
vides a better fit to the observed data. If we assume = 4.5, we get 
/* = 0.21 and 0.19 for the models with PS and ST mass functions 
respectively. 

In order to fit the observed UV luminosity function of LAEs 
of Ouchi et al. (2008) at 2 = 3.7 we need Gf = 0.1. Compar- 
ing the values of G/, we conclude that there is no strong evolution 
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in the percentage of LBGs showing up as LAEs from 2 = 3 to 
z = A. Assuming no redshift evolution in the equivalent width dis- 
tribution of Lyman-a and taking the limiting rest equivalent width 
of 45 A (as in Ouchi et al. 2008) we estimate G/ ~ 0.1 from the 
Fig. 8 of Shapley et al. (2003). However, Reddy et al. (2008) report 
an evolution in the Lyman-Q equivalent width distribution of LBGs 
between 1.9 < z < 3.4. Continuation of this trend to higher red- 
shifts will mean G/ more than 10%. Our model predictions match 
reasonably well with these observational predictions given the error 
in measurements. 

The good agreement with the data of Lyman-a LF at 2 = 
3.7 (taken from Ouchi et al. 2008) are obtained for f*ft^T = 
0.051 ±0.014 for the model with the PS mass function with best fit 
X^/dof = 0.68. The mean Lyman-a equivalent width of the LAEs 
as predicted from this model is 148 A. For the model with the ST 
mass function one needs f*ftsT = 0.050 ± 0.015 (with best fit 
X^/dof = 1.08) and the average equivalent with predicted by this 
model is 183 A. The Lyman-a rest equivalent width distribution 
predicted for the best fit model parameters is shown in Fig. 3. For 
ri = 4.5 we get /^^r = 0.25 and 0.26 for PS and ST mass func- 
tions respectively. These values are consistent with that we derived 
for 2 ~ 3. Therefore with no or minor evolutions in the physical 
conditions in the Lyman break galaxies our models reproduce the 
observed luminosity function for 3 < 2 < 4. However, from Fig. 3 
it is clear that our models predict a mild decrease in the rest equiv- 
alent width of Lyman-Q with increasing redshift. 



3.3 Luminosity functions at 2 ~ 6 

In the panels (e) and (f) of Fig. 2, we show our model prediction of 
luminosity functions at 2 '-^ 6. The observed UV luminosity func- 
tions of LBGs at 2 = 6 are taken from Bouwens et al. (2007). First, 
the required values of /./r? are 0.081±0.001 and 0.050±0.001 for 
the model with the PS and ST mass functions respectively. The cor- 
responding reduced 1-19 and 0.63 suggesting both PS and 
ST mass functions produce good fit to the data. We need G/ = 1.0 
to reproduce the UV luminosity function of LAEs at this redshift. 
Hence 100% of the LBGs are detected as LAEs at 2 '--^ 6. This 
is considerably different from 2 = 3 or 4 where only <10% of 
LBGs are detectable as LAEs. Assuming no redshift evolution in 
the equivalent width distribution of Lyman-Q and taking the lim- 
iting rest equivalent width of 25 A (as in Ouchi et al. 2008) we 
estimate G/ ~ 0.25 from the Fig. 8 of Shapley et al. (2003). This 
means that the physical properties related to the Lyman-Q emission 
have changed considerably from 2 = 3.7 to 2 — 5.7. This conclu- 
sion depends very much on the accuracy of the observed luminos- 
ity functions. While data of Shimasaku et al. (2006) is consistent 
with that of Ouchi et al (2008), there are some discrepancies in the 
fraction of Lyman break selected galaxies that are also Lyman-Q 
emitters (see Rhoads et al. 2003; Hu et al. 2004; Ajiki et al. 2003 
and Dow-Hygelund et al. 2007). We will come back to this issue in 
the discussion section. As we have been using Ouchi et al's data in 
all redshift bins we base our conclusions on their data. 

We show our model predictions for the Lyman-a luminosity 
function of LAEs at 2 = 5.7 in panel (f) of Fig. 2. The observed 
data points are taken from Ouchi et al. (2008). To fit the observed 
Lyman-Q LF one needs f*fesT = 0.044 ± 0.017 for the model 
using the PS mass function. The best fit psr degree of free- 
dom is 0.91. Therefore, even though the fraction of LAEs has in- 
creased considerably from 2 = 3 to 2 = 6, the value of ftftdi" 
in the galaxies identified as Lyman-Q emitters which characterises 
the Lyman-a escape (for fixed /*) has changed negligibly (within 
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Figure 4. upper panel : The UV luminosity functions of Lyman-o emitters 
at z = 6.5. The observed data are from Kashikawa et al., 2006. The pre- 
dicted UV luminosity functions of LBGs at 2 = 6.5 for the ST and PS halo 
mass functions are shown by solid and dashed lines respectively. We take 
the values oi f*/r] that fits the the UV luminosity functions at z; = 6. lower 
panel : The cumulative Lyman-a LF at z = 6.5. The solid Hne is for the ST 
mass function and dashed is for PS mass function. The spectroscopic (tilled 
triangles) and photometric (filled circles) data are taken from Kashikawa et 
al. (2006). The model parameters are adjusted to fit the luminosity function 
obtained from the photometric data. 



the uncertainty of the best fit values). However, this is only true 
for the model with the PS mass function. Model that uses the ST 
mass function predicts a change in the escape of the Lyman-Q pho- 
tons at < 3cr level. For this model, the best fit is obtained with 
f*fesT = 0.028 ± 0.021 (with best fit xVdof = 0.42). The 
calculated mean equivalent widths are 72 A and 75 A for model 
with the PS and ST mass functions respectively. The predicted rest 
equivalent width distribution is shown in Fig. 3. As noted above 
we see a decrease in the average equivalent width with increasing 
redshift. 



3.4 Cumulative luminosity function at 2 — 6.5 

Fan et al. (2006) have shown, based on the spectra of QSOs, that 
there is a significant increase in the IGM neutral fraction at 2 > 6. 
As Lyman-a escape also depends on the IGM opacity one expects 
a significant change in the Lyman-a luminosity function at 2 > 6. 
Kashikawa et al. (2006) have given the integrated Lyman-Q lumi- 
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nosity function and the UV luminosity function of Lyman-a emit- 
ters at 2 = 6.5. The observed data and our model predictions are 
compared in Fig. 4. We use ft/ri = 0.081 and 0.050 respectively 
for the models with the PS and ST mass functions. These are the 
best fitted values for z ~ 6 UV luminosity function of LBGs. They 
provide a good fit to the observed UV luminosity function of LAEs 
at z ~ 6 for G/ = 1. Our model predictions of the Lyman-a 
luminosity function match reasonably well with the observed data 
(bottom panel in Fig. 4). The good agreement with the data is ob- 
tained for /,/i;r = 0.054 ± 0.012 and 0.031 ± 0.015 for the 
PS and ST mass function respectively. The corresponding best fit 
per degree of freedom are 2.30 and 2.32. These two values are 
similar to those at z — 5.7. Hence, we conclude that the evolution 
in the dark matter halo mass function is sufficient to explain the 
observed evolution in the Lyman-a LF from z = 5.7 to « = 6.5 
without any major changes in other physical properties related to 
the star formation in the high redshift galaxies. 



4 CONCLUSION AND DISCUSSION 

We have built a semi-analytical model of star formation for high 
redshift galaxies which simultaneously reproduces the observed 
UV luminosity functions of LBGs and LAEs and the Lyman-a lu- 
minosity function of LAEs in the redshift range 3 < z < 6.5. We 
fit the UV luminosity functions of LBGs by changing f t,/ri while 
we adjust G /, the fraction of LBGs detected as LAEs, to match the 
UV luminosity functions of LAEs. Finally to fit the Lyman-a LFs 
of LAEs we vary /*/i^c"- The best fit values of our model param- 
eter at different redshifts allow us to probe the redshift evolution of 
properties of galaxies. In our models we make an implicit assump- 
tion that the Lyman-a emitters are a subset of a parent population 
of normal galaxies detected through Lyman break technique. 

Within the observational uncertainties, we are able to repro- 
duce the observed UV luminosity functions of Lyman-a emitters 
by simply scaling the best fitted UV luminosity functions of LBGs 
by a constant factor G/. This basically means that at a given z the 
fraction of LBGs that are seen as Lyman-a emitters is independent 
of the UV luminosity of galaxies and mass of the dark matter halos. 
Improving the errors in the UV limiinosity fimctions of Lyman-a 
emitters will allow us to investigate the possible dependence of G/ 
on the mass of the galaxies. 

The most interesting results from our study is the redshift evo- 
lution of G/. We showed that for z ~ 3.1 the well measured frac- 
tion of Lyman-a emitters among the LBGs are consistent with the 
G/ we require to fit the three luminosity function at this redshift. 
Our model fits to the observations clearly show a strong evolution 
in G/ between z < 4 and z > 5. Physically G/ at any given 
redshift will be given by the distribution in the Lyman-a escape 
among the population of LBGs. This will be governed by E(B-V), 
line of sight H I column density, velocity field in the Lyman-a emit- 
ting region and/or the duty cycle of the burst of star formation. It 
is interesting to note that even if there is absolutely no change in 
the distribution of Lyman-a equivalent width (absorption as well 
as emission) as a function of redshift one expects G/ to increase 
with z mainly because of the decrease in the liming rest equivalent 
width of Lyman-a emission in Ouchi et al's. (2008) survey. For ex- 
ample, based on Fig. 8 of Shapely et al. we expect G/ to be 0.25 at 
z ~ 6.0. From, Fan et al. (2006) we notice that the IGM transmis- 
sion decreases by at least a factor 3 between z = 3.1 and z = 5.7 
due to Gunn-Peterson optical depth. The actual change in the IGM 
optical depth in the proximity of the Lyman-a emitter is difficult 



to quantify as it depends on the ionization efficiency of the galaxy. 
Therefore, we expect G/ < 0.25 if the properties of LBGs do not 
change between z ~ 3.1 and z ^ 5.7. Thus our results giving 
G/ = 1.0atz~6 strongly support an evolution in the physical 
properties of these galaxies with redshift. 

Ouchi et al. (2008) provides luminosity functions only at 
z = 3.1, 3.7 and 5.7. From our analysis we see a sudden jump 
in G/ between z = 3.7 and z — 5.7. In order to explore whether 
this change is gradual or not, we consider few other observations 
in the intermediate redshift. At z = 4.5 Dawson et al. (2007) 
have measured Lyman-a luminosity function of LAEs. In absence 
of UV luminosity function of their sample we are unable to fol- 
low the same procedure as earlier. However, we notice that val- 
ues of Gf and f^Ji" that fit the Lyman-a LF of Ouchi's sample 
at z = 3.7 produce a good fit to the Dawson et al. data where 
as using the best fit parameters at z = 5.7 over produces the 
abundance of z = 4.5 Lyman-a emitters. There are two indepen- 
dent measurements of luminosity functions of LAEs available at 
z = 4.86 : one by Ouchi et al. (2003) and other by Shioya et al 
(2008). Ouchi et al. (2003) covers the low luminosity end of the LF 
(5 X lO"*^ < Liy„(crg s"^) < 2 X lO"*^) where as Shioya et al 
(2008) covers the high end (8x10*^ < Liy„(crgs"^) < 4x10*^) 
with slight overlap between them. The Ouchi et al. (2003) mea- 
surements are consistent with Gf = 0.1. However, if we also con- 
sider Shioya et al (2008) data, Gf could be as large as 0.3. Note 
that the completeness of the sample is always an issue in this case. 
Hence more observations are needed in this redshift range in order 
to probe in detail how G/ increases to unity by z = 5.7. 

Unlike at z ~ 3.1, the luminosity function of Lyman-a emit- 
ters obtained by different groups for z ~ 5.7 disagree up to a factor 
5 (see Rhoads et al. 2003; Hu et al 2004; Ajiki et al. 2003; Mu- 
rayama et al. 2007; Shimasaku et al. 2006; and Ouchi et al. 2008). 
The difference could be due to differences in the colour selection 
criteria used in the narrow band survey and the depth of the broad 
band photometry. Dow-Hygelund et al. (2007) have found that 
only 30% of the Lyman break galaxies at z ~ 6 selected through 
i-dropout selection show Lyman-a emission with rest equivalent 
width > 20 A. It is also important to remember that while the nar- 
row band imaging picks object within very narrow redshift range 
the broad band colour techniques pick objects over a much wider 
redshift range. Incompleteness levels in these two types of surveys 
are also very different. Dow-Hygelund et al. (2007) have shown 
that the i-band selection misses considerable nimiber of Lyman-a 
emitters at z < 5.8. On the other hand, the narrow band technique 
of Ouchi et al (2008) picks object at z = 5.70 ± 0.01. After tak- 
ing into account this effect Dow-Hygelund et al (2007) conclude 
that up to 40% of the i-dropout galaxies could be Lyman-a emit- 
ters. From our models we find the redshift evolution between the 
mean z of LBGs and LAEs will account for an additional 10% in- 
crease in Gf. Even after taking into account all these effects one 
needs G / to be factor 2 higher to explain the available observed lu- 
minosity functions. Thus we can conclude that there is an increase 
in G/ as a function of z but to get the actual amount we need lot 
more observations at z > 5. Recent results from the narrow band 
survey of Lyman-a emitters at z ~ 4.7 by Shioya et al (2009) are 
also consistent with increasing value of G/ with increasing z. As 
Lyman-a escape depends on the amount of dust and gas kinemat- 
ics, the higher value of G/ implies that on an average the ISM of 
z > 5 galaxies are less dusty, more clumpy and having complex 
velocity field making the escape of Lyman-a photons easier. 

Further, the evolution in the observed Lyman-a LF at z > 5.7 
can be imderstood as evolution in the number density of the dark 
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matter halos arising from the structure formation model with mod- 
est change in the physical properties of these galaxies. This is inde- 
pendent of the form of the halo mass function we assume. Dijkstra 
et al. (2007) have arrived at same conclusion in the evolution of 
luminosity functions for 2; > 5.7 while considering no evolution in 
the IGM transitivity in this redshift range. 

Our best fit models at different redshifts show that average 
Lyman-a equivalent width decreases with increasing redshift. This 
is contrary to some preliminary observational results that suggest 
an increase of equivalent width with increasing redshift (Grove et 
al. 2009). This result needs to be confirmed with larger number of 
spectroscopic data. In our model it is possible to get such a trend by 
allowing the initial stellar mass function to vary with redshift (see 
Fig. 1). Also, relaxing our assumption that the equivalent width is 
independent of galaxy mass will have some effect on the equiva- 
lent width distribution. Indeed such mass dependence of equiva- 
lent width distribution is indicated by observations of Ando et al. 
(2006). 

There are a number of other attempts to fit the U V and Lyman- 
a luminosity functions using galaxy formation models. Kobayashi 
et al. (2007) using their hierarchical galaxy formation models fitted 
the luminosity function of Lyman-a emitters by varying the es- 
cape fraction. However, according to their models all LBGs would 
be detected as Lyman-a emitters. Mao et al (2008) fitted luminos- 
ity functions using semi-analytic models that compute E(B-V) and 
relate it to the escape fraction of Lyman-a photons. In this model 
preferably low metallicity dust free galaxies will be seen as Lyman- 
a emitters. However, recent observations suggest that the Lyman-a 
emitters need not be confined to primordial low dust populations 
(Pentericci et al. 2008; also see Scannapieco et al. 2003, Fynbo et 
al. 2003, Dawson et al. 2007). Nagamine et al (2008) used the hi- 
erarchical structure formation models to fit the Lyman-a emitters 
assuming a normal galaxy is a Lyman-a emitter for a brief period 
of time (duty cycle argument). They find the duty cycle increases 
with increasing redshift as we find for G/. 

It is important to realize that high redshift luminosity functions 
are based on deep field observations covering small volumes. The 
effect of cosmic variance may be large. The UV luminosity func- 
tions used here for LBGs are mainly based on photometric data 
with large redshift uncertainty. Therefore, more observations are 
needed to get a clearer picture on the evolution of physical proper- 
ties of the galaxies. In the case of modelling, one requires a clear 
physical model for G/. It is possible that simple ideas of duty cy- 
cle based on dust properties may not be sufficient since the veloc- 
ity field in the Lyman-a emitting regions may play an important 
role. Indeed, all the high- 2: LBGs show signatures of outflows that 
can enable easy transport of Lyman-a photons. Thus, physical un- 
derstanding of G f based on a dynamical model (e.g Verhaimne et 
al. 2008) that will also fit the liuninosity functions is the next step 
in this subject. Such models may also explain the observed wide 
spread in the rest equivalent width distribution. 
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